Stability of tardyons and tachyons in the rotating and expanding Universe 
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In the present paper we analyze the spectrum of quasinormal modes for massive scalar and Dirac 
fields, allowing for both tardyonic (/i 2 > 0) and tachyonic (/i 2 < 0) masses, in the expanding and 
rotating cosmological background. The spectrum found shows a number of peculiar features, which 
are absent in the Minkowski space-time. A hypothetical particle that moves faster than light, a 
tachyon, is known to be classically unstable in the Minkowski space-time. This instability has its 
analog at the quantum level: small vacuum fluctuations of the field lead to the unbounded growth 
of the amplitude, so that the appearance of the real tachyons in the spectrum means that there is 
catastrophic instability in the theory. It has been conjectured a long time ago that possibly the 
lightest particles with nonzero mass, the neutrino, may be a tachyon. Here we shall show that in 
the rotating and expanding Universe tachyons are stable if their mass is less than some constant, 
which is related to the Universe's rotation and expansion scales. The current upper bound on the 
rotation scale gives us a very small upper bound on the tachyon's mass which is many orders less 
than the mass of the electron. This might explain why only very light particles have the chance 
of becoming tachyons. It is shown that the spectrum of the "normal" (fi 2 > 0) Dirac field has a 
discontinuity as a function of the z-component of the wave vector k z at k z = 0. 
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I. INTRODUCTION 

The spectrum of classical fields in various cosmological 
space-times is a well-established area, which is a starting 
point for quantization of the fields in the curved back- 
grounds. By now it includes a detailed knowledge of spec- 
tra for massless scalar and Dirac fields [3-0 i n the rotat- 
ing (but not expanding) Godel-like universes [1, H| ■ The 
massless scalar field in the expanding and rotating back- 
ground was considered by Panov [6j. Yet, all these works 
were devoted solely to "normal" nontachyonic fields. 

The first comprehensive study of would be particles 
moving faster than light, i.e. tachyons, was done as 
early as 1962 Q. Since that time there have been a 
number of other theoretical studies on the various prop- 
erties of these hypothetical particles in quantum field 
theory, string theory, and gravitation Q. One of the 
main problems with tachyons is their instability. In- 
deed, one can easily see that in the Minkowski space-time 
ds 2 = dt — dx 2 — dy 2 — dz 2 , the Klein-Gordon equation 
with the negative values of fj, 2 



+ ^ = (1) 



allows for the normalizable exponentially growing modes 
(f> (x e~ lujt , where ui = Re(u)) + ilm(uj) : Im{uj) > 0. 
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This means that vacuum fluctuations of the field will 
grow without the bound leading to the instability and 
decay of the tachyonic state. Strictly speaking, the above 
instability creates unsurpassable difficulties in quantiza- 
tion, so that the "tachyon" has never been considered se- 
riously by the majority of researchers as a real quantum 
particle, but rather like a synonym of the onset of in- 
stability in the system. Tachyons also play an important 
role in cosmology, where they are used in the spontaneous 
symmetry breaking 

Theoretically, the tachyonic nature of the neutrino was 
proposed in the 1980s [10], and once one accepts the 
conception of the tachyon, it would be natural to ex- 
pect a number of different tachyonic particles with vari- 
ous masses. Yet, today we observe only single particles, 
neutrino, which potentially could be a tachyon. Though, 
there are also arguments against the tachyonic nature of 
the neutrino, to which we shall possibly contribute one 
more argument here. If one supposes the tachyonic na- 
ture of a superluminal neutrino, then a straightforward 
question is why we observe only a single candidate for 
the tachyon, with a mass that is apparently much smaller 
than the electron mass and why "heavy" tachyons cannot 
be observed? In this work we shall show that only very 
light tachyons can be classically stable in a number of 
rotating cosmological backgrounds. At first glance, this 
could answer both questions posed above: the instabil- 
ity and the selectivity in favor of the lightest particles. 
The Universe is supposed to be expanding, so that the 
initial rotation should have mostly been damped today. 
Therefore, as only a very slow rotation of the Universe 
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is allowed by current observations, the stability is possi- 
ble only for those tachyons that are many orders lighter 
than the electron. Yet, we observed that the threshold of 
instability occurs at a tiny mass which seems to be too 
small to provide superluminality of the Standard Model 
particles. 

Previously, one of us reported on the stability gap 
for tachyons in the pure Godel-like nonexpanding space- 
times Here we shall generalize this result to a more 
realistic situation of expanding, rotating Universe and 
shall make estimates for the upper bound of tachyon's 
mass from the current upper bound on the rotation rate 
of the Universe and the Hubble constant. 

Another aim of our research, which is not connected 
to supeluminality issues, is analysis of spectra of "nor- 
mal" (tardyon) fields in the expanding and rotating Uni- 
verse. We have found a number of distinctions from the 
previously studied cases of scalar fields in these back- 
grounds, which may raise new problems under quantiza- 
tion of these fields. 

The work is organized as follows. Sec. [H] is devoted to 
the stability of a scalar tachyonic field, and Sec. Mil deals 
with massive Dirac tachyons. Sec. IIVI summarizes the 
spectrum of tardyonic scalar and Dirac fields. In Sec. |V] 
we discuss the obtained results and open problems. 



II. SCALAR FIELD 

A wide class of rotating and expanding cosmological 
solutions can be described by the Godel-like metric 



da' 



(dt + a(t)e ax dy) 2 -a(t) 2 (dx 2 



B 2 e 2ax dy 2 



dz 2 ) , 
(2) 

which includes both causal and noncausal generalizations 
SEMI of the classical Godel solution j|. The Godel 
metric is reproduced if B 2 = 1/2 and a(t) = 1. 

Although here we shall use the Cartesian coordinates, 
let us note that the metric ^ can be rewritten in the 
cylindrical coordinates (t',r,cf>,z) as follows: 



ds 2 = (^a(i')sinh 2 
a 2 (t') ( dr 2 



dz 1 



sinh (ar) 



(3) 



If the following inequality is fulfilled 

B 2 -1 



1 



B' 2 



■ sinh' 



(?) 



>0, 



or simply, 



B > 1, 



(4) 



(5) 



no closed time-like curves are allowed for the space-time. 
This is an example of the causal cosmological solutions 
suggested in One can see that value of B does not 
change our estimations for the order of the tachyon mass, 



so that, qualitatively, our conclusions will not depend on 
the causality of the background space-time. 

Following Q we introduce the new coordinates r and 
•q instead of t and x, 



dr = 



dt 

W)' 



V 



In terms of these new coordinates the metric takes form, 
dy" 



ds = a(r) 



dr- 



2 drf B 2 dy 2 2 
dz 



<J 



B 2 a s ( T ) 
a 2 rf 



(6) 



We can separate the variables in the Klein-Gordon equa- 
tion (P), using the following ansatz 

cj ) (r,r 1 ,y 7 z) = e^y + ^ z R(T,f ] ). (7) 

Then, the wave-like equation can be written as follows: 



d 2 R l-B 2 (d 2 R a'{T)dR\ 2ik y dR 



drf B 2 a 2 r) 2 \ dr 2 a(r) dr J B 2 a 2 rj dr 
( k 2 d | k 2 z +a 2 (r)n 2 \ R a'(r) 2ik y R () 



(8) 



B 2 



7] 



a 2 a(r) B 2 a 2 r\ 



Let us be limited by an adiabatic expansion of the uni- 
verse, i.e. we shall consider the radius of the universe 
a(r) and the Hubble parameter H as constants: 



1 da a'(r) 
a[t) at a (j) 



(9) 



Then, following [6J , we introduce the frequencies of the 
proper oscillations of the universe as 



R(T,r,)=^2e- i ^N u (ri), 



(10) 



where N u (j]) satisfies the following equation 



2k y w + 2ik y aH k 2 + a 2 p? 



B 2 



B 2r q 



(11) 



(1 - B 2 ){lo 2 + iujgH) \ NUv) 
B 2 V 2 J a 2 

Eq. (|11[) has a regular singular point at r\ = and an 
irregular singular point at r\ = 00 with the asymptotical 
behavior 



N u cx exp ( ± 



k v r\ 



Let us notice, that the asymptotically Godel space- 
time, as well as the asymptotically anti-de Sitter (AdS) 
one, is not globally hyperbolic, so that mathematically 
strict analytical proof of the stability should also include 
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the detailed consideration of the well-posedness of the 
initial value problem. We shall study only normalizable 
solutions, i.e. 



-g drj < oo 



^f-dn < oo. (12) 



In order to satisfy convergence conditions at both end- 
points of the integral we require 



N u (rj — > oo) oc exp 



p> 1/2 

_\ky[n 

Ba 



In terms of the new variable £, 



\ky\T] 



Ba ' 

the equation takes the following form 



where 



J 
K 



m(uj + iaH) 
Ba ' 

fc| + a 2 pi 2 (I- B 2 )(w 2 +iLoaH) 
a 2 a 2 B 2 



(13) 
(14) 



Q + | + f)^(6 = o, (is) 



(16) 
(17) 



and m = sgn{k y ). Taking into account the asymptoti- 
cal behavior of the wave equation, we introduce a new 
function P(£), 



(18) 



Choosing q such that K = q{q — 1), we find that the 
function P(£) satisfies the Kummer's equation 

£P"(0 + (2q $P'(0 = (q + J)P(£). (19) 

The only solution, which satisfies (fT3")) as £ — > 0, is 

_ f M(g + J, 2o, £), if #e ( g ) > 1/2; 

UJ ~ 1 ^- 2 "M(J - q + 1, 2 - 2q, £), if i?e ( 9 ) < 1/2; 

(20) 

where M(a,/3,£) is the confluent hypergeometric func- 
tion. 

As £ — > oo, M(a, /3, £) grows exponentially and violates 
(|14l) unless a is a nonpositive integer. Thus, 

I q + J, if Re (q)> 1/2; 
JV - \ J-q + 1, if Re (q)< 1/2. ^ 

where N = 0, 1, 2, 3 . . . Finally, we find 

if= (JV + J)(N + J+l). (22) 

A general solution to (|22[) gives us an analytic expression 
for the frequency to: 



(l + 2N)a TT (l + B 2 ) a 2 + 4(k 2 +(i 2 )-a 2 H 2 (B- 2 -l) 2 fiaH(l + 2N)a (l + 2N) 2 a 2 \ ,„ „ 
^5 mi? ^B^ ± V " 4 + { 2B^ + 4 ) {B l) - 

(23) 



When B 2 = 1/2, H = and a = \/2fi, we have the 
nonexpanding Godel space-time. When the expression 
under the square root is positive in the above equation, 
oj is pure real, so that the tachyon is stable. Thus, the 
upper bound on the tachyon's mass reads 

M < n, (24) 

which coincides with Eq. (13) of (Tlj . 

When B = 1, the expression for uj is formally reduced 
to the one in KjJ: 



(1 + 2N)a ± y/a 2 + 4(k 2 + a 2 y?) 

uj — — 

2m 

There are unstable modes in the spectrum if 
a 2 +4a 2 f i 2 > -4a 2 i? 2 , 



allowing for tachyons with very small mass 

where Qh = Q/H = a/\2aH is the rotation parameter 
measured in units of the Hubble constant. Estimations 
in [15[ show that can be equal to 74. 

The Hubble constant is of order 10 -10 i/ear -1 , giving 
an estimation for the tachyon mass 

N £ W^year^h/c 2 - lO^kg ~ 10" 36 m e , (26) 

where m e is the electron mass. 

The estimation for the upper limit on the Universe's 
rotation certainly crucially depends on the cosmological 
model under consideration (see for instance a recent re- 
view on the observational aspects of the global rotation 



4 



[III). Therefore, the upper limit for the global rotation 
that we have used is most physically relevant as it deals 
with the same cosmological model. 

Apparently, the most authoritative approach to finding 
an upper constraint on the Universe's rotation is based 
on the Cosmic Microwave Background Anisotropics 
(CMBA) data of NASA, owing to their precision. 
Constraints for nonflat cosmological models, including 
the Godel-Friedman ones, are in the r ang e 10~ 9 
10~ 8 rad/year [l5l - [T8| . In particular, in [l8l ] the anal- 
ysis of the so-called Sachs- Wolf effect (the gravitational 
redshift of photons in the cosmic microwave background, 
which makes the cosmic microwave background spec- 
trum uneven) leads to the 10~ 9 rad/year constraint for 
the global rotation. If one considers this constraint of 
10 rad/year as physically relevant, then the estimated 
upper limit on mass stability will be one order less than 
in Eq. 1[2B]>. 



and ef- is the tetrad field which satisfies 



til/ ab n v ilv ab a v 

9 = V &a e b > 9 = V e a e-b ■ 



One can choose the tetrad field as follows 



Ba(r) 



/ B 

Bar] 

-1 i) 

V JS 



(28) 



Here we shall use the standard representation for the 
Dirac matrices 



(29) 



i \ 




( 


ia 1 \ 


-i I ' 




V ia 1 


) 



III. DIRAC FIELD 

The Dirac equation in a curved space-time has the form 
-i 7 a e/($M + r M $) + ^ = 0, (27) 
where j a are Dirac matrices 

{ 7 V} = V, 



where m are standard Pauli matrices (i — 1,2,3). 
The spin connections T M are defined as 



1 



■b ;n' 



(30) 



For the space-time under consideration it takes the fol- 
lowing form 



hW] _ ^M [7 °, 7 *]) ; r, - Jg ( b Wi + ^ (b [7 », 7 .] + b Wl)) , 



8r] 

As for the scalar field, we shall consider the adiabatic approximation for the Universe's expansion and assume the 
following ansatz: 

$(r, 77, y, z) = e-^r+ikyv+iKz^^^ ^ 

Replacing $ in the Dirac wave equation by pip , and after some calculations, one can reduce the wave equations to 
the following the matrix form: 

*M = ( M Q ) whcrc M = ^ v ( a, °0 ) ' ( 32 ) 



P 



Q 



2(Ba ~ 2k y T] - 3iaH - 2cu) ~6BaH - AiBk z + AiBoj + i 

ABari \ ~6BaH + AiBk z + AiBu: - ia 2{Ba + 2k y r] + 3iaH + 2uj) 

2(Ba - 2k y i] - 3iaH - 2u) 6BaH - AiBk z - AiBuj + ia 



ABa-q \ 6BaH + AiBk z - AiBuj + ia 2(Ba + 2k y ij + 3iaH + 2uj) 



Choosing the 4-component spinor as 



m = U + -7 3 



Mv)] , k z =k z Ul + ^--l\ (fc z 2 =aV, if k z = 0), (33) 



5 



we are able to reduce the matrix equation to the second-order differential equations for each component 



k 2 y 2k y u ± k y Ba AF \k z ) 2 -a 2 (1 - B 2 )C? \ i\)(r\) 



B 2 



B 2 r) 



4?7 2 



B 2 T] 2 



0, 



(34) 



where Co = u> + 3iaH/2. For left-handed components 



a 2 /i 2 + k 2 



a 2 n 2 + k 2 4B 



lim F{k z ) = k z 



AB J 



In the limit H — > 0, \i — > we reproduce formula (102) of [l[ (one should take B 2 = a = 1/2). 
Again, we study only the normalizable solutions 



J *7o 



<&yj—g drj < oo 



— drj < oo, 



(35) 



(36) 



which imply the same requirements for the endpoints as for the scalar field (|13|) . (|T4l 
Eq. (|34[) has the same form as (fl"9)) with the coefficients J and K given by 



muj 1 

J = ± -, 

Ba 2' 



K = 



F 2 {k z ) 1 (1-B 2 )Cj 2 



a 2 B 2 



From (l22l) one finds 



(37) 



(2AT+l±l)a ^(27V + 1 ± l) 2 (B- 2 - l)a 2 + AF 2 {k z 
U) — — — zb 



3iaH 



2Bm 



(38) 



The "-" sign for N = violates the conditions of convergence of the Dirac field norm (see [l[ for details), so that we 
are able to write both solutions as 



± y/fN + l) 2 (B- 2 - IW + F 2 (k z ) - 
Bra I 



(39) 



r 



For B = 1 the spectrum of the wave-like equation takes on the rotation parameter. The instability condition be- 
comes 



the following simple form 
(N+l)a 



±F(k z 



3iaH 



771 2 

Then, for B = 1 the instability condition is 

, 3aH 
Im{F(k z ))>—. 



Since the imaginary part of F(k z ) reaches its maximum 
for k, = and 



Im y/(N + 1) 2 (B- 2 - l)a 2 + F 2 (k z ) ) > 



3aH 



(41) 



From the formula 



\Im (Vx) 



x\ — Re (x) 



(Re (x) > 0), 



max(Im (F(k z ))) = Im (F(0)) = Im (afi) , we conclude that N = corresponds to the most unsta- 

ble mode, for which we find (for /i 2 < 0) the following 
the bound for fermion mass does not depend on the Uni- expression: 
verse's rotation in this case, and, 

2.H 



< 



(40) 



a 2 (17- 16S 2 ) 



+ a 2 /i 2 



However, for B 2 ^ 1, and, in particular, for the Godel 
case (B 2 = 1/2), the bound on the tachyon mass depends 



16S 2 

q 2 (17- 16B 2 ) 
16B 2 



a 2 a 2 [i 2 



+ a 2 fi 2 I > 



4B 2 

9a 2 H 2 



6 



Finally we find that the tachyon is stable if 



< 




(42) 



i8B 2 /n 2 H ' 



For the Godel-like solution (B 2 = 1/2) and the quick 
rotation scenario (fin 3> 1), we obtain the tachyon mass' 
bound of half order higher than the Hubble parameter 



I Ml 



< 



9H 



(43) 



Let us notice, that if B 2 > 1 the expression inside 
the square root in (|23|) . ([38} can always be negative for 
sufficiently large TV, even for the massless case. Thus, for 
B 2 > 1 tachyons (both scalar and Dirac) are unstable, 
yet the instability occurs for the normal matter as well, 
making this case quite exotic. In summary, noncausal 



space-times with B 2 < 1 and causal with B 2 = 1 allow 
for the tachyon's stability, while causal space-time with 
B 2 > 1 is a quite exotic one, because not only tachyonic 
matter, but also the normal one is unstable. 



IV. QUASINORMAL MODES OF 
NON- TACHYONIC FIELDS 



First let us note, that solutions of (|2"21 are not always 
also the solutions of (|2"Tj) . Because of the requirement for 
the real part of q (HH), one can see that J must always 
have a negative real part, which implies that the real 
part of u) has an opposite to the m sign and always has a 
minimal value. That is why if one studies both the real 
and imaginary parts of w, he must discard solutions that 
do not satisfy pTjl . 

When [i 2 > we can prove that Eq. (|2"Tj) is not satis- 
fied for all of the modes of (j2"3")l for which the sign of the 
real part of the square root is — m. Indeed, [24j, 



Re {mu) > _M + Re >^<r-" + (HS^Da + (£j -2 — i) I m 



2B \ \j 4B 2 \ 2Bm 

(l + 2/v)o , Re gi + 2N)a + iaH\ = Jl + 2N)a (1 _ > Jl + 2N)a 



2B \ 2 2Bm) 2 \B 

from which it follows that Re (J) > —N— |. 

Therefore, for fi 2 > the quasinormal spectrum is 



(l + 27V)a (1 + B 2 ) 1 a 2 + 4(k 2 +ii 2 )-a 2 H 2 (B- 2 -l) 2 fiaH(l + 2N)a (l + 2N) 2 a 2 \ 

U = TT-^ ™H „ - —J + —— + (B --1) 



2mB ' " 2B 2 m\j 4 V 2Bm 

(45) 

Similarly, we can prove that for the Dirac field 



(JV + l)q ^(N + l) 2 (B~ 2 -l)a 2 +F 2 (k z ) 3zaH ^ 2>Q 



Bm m 
This was also shown for a particular value of B, B 2 = |, in [l[ 



One can see, that for H ^ the frequencies have an 
imaginary part of order aH, which corresponds to the 
decay of the amplitude due to expansion of the Universe. 
The real part of the frequencies form qualitatively the 
same picture as that for H = (see Figs. I1I2[) . 

From Fig. [2] we observe that quasinormal spectra are 
not the same for particles of opposite chiralities. The 
spectrum for the right-handed chirality can be obtained 
by the replacing k z — > —k z in (l3"5j) . Since for massless 
particles the helicity is an invariant quantity, we observe 
a continuous dependence of the QNMs on k z (left panel 
of Fig. [2]). Discontinuity for the massive particles of the 
same chirality as k z — > (right panel of Fig. [2]) reflects 



a nontrivial dependence of QNMs on helicity, which is 
not an invariant quantity in this case. The left and right 
limits as k z — > correspond to particles with different 
dominant helicities that are described by the same left- 
handed chirality. If one considers opposite chiralities for 
the positive and negative values of k z , implying domina- 
tion of the same helicity, the figure is symmetric for both 
massive and massless particles and has no discontinuity 
when /i 2 > 0. 



Note, that normalizable solutions exist only when k y 7^ 
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FIG. 1: Left panel: real parts of Dirac quasinormal modes as functions of F(k z ) for B = 1/2. Upper (positive real part) curves 
correspond to m = — 1, lower (negative real part) correspond to m = +1. Right panel shows k z as a function of /i for which 
F(k z ) is zero. 




0. For k y = 0, Eq. (JTTJ) reads 

K(V) = ^N U ( V ), (47) 

whose general solution is 

N u (r,) = ^n(C + rf + C-r,- p ), 

where p 2 = K + \ . 

The only nondivergent solutions is constant when p = 
±5 and K = However, for these solutions the 
integrals (fl"2|) and (J36J) diverge. 



V. CONCLUSIONS 

In this work we have shown how the cosmological rota- 
tion and expansion can influence the stability of tachyons 
and the spectrum of classical tardyon fields. The stability 
gap for tachyons was expected from the similarities be- 
tween the spectrum of asymptotically Godel [l9l - l2l1 | and 



AdS space-times. Indeed, in (2^, [2l[ it was noted that 
the spectrum of black holes in the 5-dimensional asymp- 
totically Godel space-time, is qualitatively similar to the 
normal modes of pure Godel space-time and is quite dif- 
ferent from the quasinormal modes of the black hole in 
asymptotically flat space-time. In the limit of the vanish- 
ing black hole radius, quasinormal modes approach the 
normal modes of the pure Godel space. The same behav- 
ior takes place for black holes in AdS space-time [22|, [23[ , 
where the anti-de Sitter radius R plays the role of the 
Universe's rotation rate ft. 

Here we have found that the spectrum and stability 
of Dirac tachyons is quite different from the scalar one: 
while the scalar tachyons do not need expansion for the 
stability gap, that is, the gap exists for rotating and not 
expanding models, Dirac tachyons are unstable in the 
nonexpanding cosmological models. Yet, for both scalars 
and neutrino tachyons the coupling of rotation and ex- 
pansion can enlarge the stability gap, making thereby 
cosmological rotation an important stabilizing factor. 

Finally, since we have shown that the upper bound on 
the tachyon mass in all of the above cases is quite small, 
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and is definitely many orders less than the electron's 
mass, one could suppose that this could explan why today 
only the lightest particles (though with nonzero mass, 
such as neutrinos) could be tachyonic. 

An important issue which was beyond our consider- 
ation is quantization of the considered fields in the ex- 
panding Godel-like backgrounds. One can show that the 
spectrum of bound states considered here do not form a 
complete set, so that the quantization of massive fields in 
such space-times is a nontrivial question, not to mention 
the tachyonic nature of the fields. Yet, if such a quanti- 
zation is performed, the effort will pay off, allowing for a 
quantum description of stable tachyons. 

In addition, we have described the spectra of tardyon 
fields in detail. The spectrum of the massive scalar or 
Dirac fields do not form a complete set and do not have 



unstable normalizable modes. Frequencies for a massive 
Dirac field (as well as a massless one) are not symmetric 
relative to the k z axis, where k z is the component of the 
wave vector. A new feature which that we observe for the 
massive Dirac field is the discontinuity of the frequency 
ms a function of k z at k z = 0. 
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